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A method to create paired atom laser beams from a metastable helium atom laser via four- 
wave mixing is demonstrated. Radio frequency outcoupling is used to extract atoms from a Bose 
Einstein condensate near the center of the condensate and initiate scattering between trapped and 
untrapped atoms. The unequal strengths of the interactions for different internal states allows an 
energy-momentum resonance which leads to the creation of pairs of atoms scattered from the zero- 
velocity condensate. The resulting scattered beams are well separated from the main atom laser in 
the 2-dimensional transverse atom laser profile. Numerical simulations of the system are in good 
agreement with the observed atom laser spatial profiles, and indicate that the scattered beams are 
generated by a four- wave mixing process, suggesting that the beams are correlated. 
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Sources of matter waves gained a dramatic improve- 
ment with the achievement of Bose-Einstein condensa- 
tion (BEC) in dilute gases and the development of the 
atom laser [TJ [2] . Like optical lasers before them, atom 
lasers can produce Heisenberg-limited beam profiles [3] 0] 
and promise high spectral density through their dramat- 
ically lower linewidth [5 . Another exciting possibility 
resulting from having such a coherent source of atoms 
is the generation of non-classical matter waves through 
entangled beams. Such entangled beams are useful for 
tests of quantum mechanics, and are required to perform 
Heisenberg-limited interferometry [BJ[7]. In this Letter, 
we show that the asymmetric scattering rates between 
internal states of metastable helium (He*) cause well- 
defined peaks in the output of an atom laser. These peaks 
are due to a four- wave mixing (FWM) process, and are 
experimentally demonstrated. 

A nonlinear process is required to produce entangle- 
ment, and one of the advantages of atomic systems over 
optical systems is that there are strong inherent non- 
linearities due to the atomic interactions, although these 
interactions can also lead to complications. These nonlin- 
earities allow certain analogs of nonlinear optical experi- 
ments such as FWM and Kerr squeezing to be performed 
directly in the atomic sample [8|. Both of these produce 
entanglement in optical systems. FWM in a trapped 
BEC has been demonstrated experimentally in configu- 
rations where three distinct momentum states generated 
a fourth [9] , and where two momentum states generated a 
pair of correlated atomic beams [10]. These experiments 
demonstrated that the output was phase coherent, but 
the correlation properties were not measured. More re- 
cently the pair correlations in a spontaneous scattering of 
two colliding condensates were measured using the single 



atom detectors available for He* atoms [TT] . 

Using these existing sources of entangled pairs of atoms 
for interferometric experiments will be complicated by 
the high densities of the sources, where the nonlinear- 
ities that generated the correlations ultimately degrade 
the long term coherence of the sample. While recent 
experiments have increased the coherence of atom in- 
terferometers by several orders of magnitude by reduc- 
ing the nonlinearities with a Feshbach resonance p~2j [13] , 
this precludes the production of entangled pairs. In our 
scheme the nonlinear interactions are used to drive FWM 
in the magnetically trapped condensate, but the result- 
ing untrapped beams that propagate in free space are di- 
lute, potentially avoiding the decoherence problem. Us- 
ing atoms in the untrapped state also makes the beams 
insensitive to magnetic field inhomogeneities. We show 
that pairs of beams can be produced simply by the pro- 
cess of radio frequency (RF) outcoupling from a He* 
BEC, without the need for Feshbach resonances, optical 
traps or scattering pulses. Unlike the previous methods, 
which required pairs of atoms travelling at high kinetic 
energies as a source, this process involves scattering be- 
tween atoms initally in the same zero momentum state 
to create states with non-zero momentum. The energy- 
momentum resonance comes from the mean field con- 
ditions that are obtained during outcoupling from the 
condensate. Semiclassical and field theoretic simulations 
of the experiment show that the beams are generated by 
the same FWM process that generated entangled atom 
pairs in the earlier experiments. 

Our experimental setup for creating a He* BEC has 
been reported elsewhere [13] . Briefly, we use a cryogenic 
beamline to produce a He* magneto-optic trap (MOT). 
Atoms are extracted from this low vacuum MOT into 






a high vacuum MOT via an LVIS + setup. Atoms are 
transferred from the MOT into a BiQUIC magnetic trap 
(with axial trapping frequency of 55 Hz and radial trap- 
ping frequency of 1020 Hz), where a BEC is produced 
using forced RF evaporative cooling. Using this setup 
we are able to create almost pure BECs containing up to 
5 x 10 6 atoms. To produce our atom laser beam, we use 
RF photons to spin flip the BEC atoms from the m = 1 
magnetically trapped state to the m = untrapped state. 

After outcoupling, atoms in the atom laser beam fall 
under gravity for a distance of 4 cm until they strike a 
double stacked multi-channel plate (MCP). We image 
the phosphor screen with a CCD camera with a reso- 
lution of approximately 150 /im at the MCP. To remove 
any nonuniformities caused by spatial variations in the 
gain of our MCP, we divide all our images by a flat field 
image, produced by dropping atoms from a MOT onto 
our detector. Since the MOT temperature is of order 
~ 1 mK, the spatial profile of the MOT uniformly illumi- 
nates our MCP. Although m = — 1 atoms are produced 
by the outcoupling, especially for high RF powers, they 
are in general accelerated away from our detector by the 
magnetic trap field. Those that are accelerated towards 
the detector do not show up in our images since they 
arrive much earlier than the CCD trigger time. 

Figure [T] shows the dramatic change in the atom laser 
spatial profile when the condensate undergoes the reso- 
nant FWM process. In the case of low outcoupling Rabi 
frequencies (upper row) the FWM process does not occur 
and we see the usual double peaked He* atom laser pro- 
file [15] . When the Rabi frequencies are high enough to 
initiate the FWM process (middle row) the atoms in the 
atom laser beam are scattered to form a halo around the 
main atom laser beam. Due to conservation of energy 
the outer diameter of this halo corresponds to a maxi- 
mum kinetic energy given by the chemical potential. As 
well as the ring structure, four peaks are observed on the 
outskirts of the profile. Each peak contains 11-14% of 
the total out coupled atoms, in agreement with the value 
of 11% obtained from the simulation results shown in the 
lower row of Fig. [I] 

These peaks arise from two momentum-correlated 
cones of particles scattering out of the condensate and 
falling to the detector, as shown in Fig. [2] The cones 
themselves are generated via a FWM process that pop- 
ulates momentum modes lying along the weak trapping 
axis, which expand to rings due to the mean field re- 
pulsion in the tight trapping direction as they leave the 
BEC. Atoms in these cones then fall under gravity onto 
the detector, where the time-integrated flux is converted 
to a spatial density distribution, and each momentum 
ring appears as a double peak. The background halo is 
produced during the initial switch-on of the atom laser 
during which the peaks produced by the FWM process 
sweep in position from the main atom laser profile to- 
wards their final position in Fig. [I] 
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FIG. 1: (color online) First two rows show experimental atom 
laser spatial profiles observed on our MCP 4cm below the 
trap, in a 3-D rendering (left) and the 2-D image (right). 
Both sets of data were taken for an outcoupling detuning of 
2 kHz, however the Rabi frequency is increased by an order of 
magnitude between the two sets. The upper row shows the 
usual He* atom laser [15] , while the middle row demonstrates 
the appearance of the resonant scattering peaks. The bottom 
row is the result of a simulation of the second experiment. 
The weak axis of the trap is aligned along the vertical axis of 
the right hand images. 



In order to understand how these peaks are produced 
we carried out a detailed numerical simulation of the ex- 
periment. The master equation we used to describe the 
experiment is given by 
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FIG. 2: (color online) Schematic diagram of experimental 
configuration. The four- wave mixing process generates mo- 
mentum components along the weak trapping axis, which get 
outcoupled to two cones. After out coupling from the conden- 
sate, atoms fall 4 cm to the MCP detector. 



where 

H=Y.jf *}(x) ("^ + j Wx)) *j (x) dx 

+ Eij I *t( x )v^(*ij+i + ^-i)M^(x) dx (2) 
+ Eij U ~f I *l(x)*](x)^(x)*,(x) dx, 

and where T>[c]p = cpc^ — | {c^cp + pc^c) , \£ m are 
the field operators for the magnetic sublevels ra, fi = 
350 Hz is the Rabi frequency of the RF outcoupling, 
Vtrap(x) = t>M (^x 2 + c^ 2 ?/ 2 + oo 2 z z 2 ) is the trapping po- 
tential, Uij = 4:7rh 2 aij /M is the nonlinear interaction 
strength, dij is the s-wave scattering length between the 
magnetic sublevels i and j, and iC^ po1 ^ is the Penning 
ionisation rate for an unpolarised thermal sample of He* 
[16] . The scattering lengths for He* are a\-\ = 3.33 nm, 
a o,o = 5.56 nm and 7.51 nm for all other combinations of 
internal states [17] • The Penning ionisation loss term in 
Eq. ([!]) was obtained by assuming a position-independent 
loss for the molecular state with total angular momentum 

of 0, = ^(2* 1 (x)*_ 1 (x)-*o(x)*o(x)) ) 

which is the dominant contribution to Penning ionisation 
at condensate temperatures [16]. The rate constant was 
determined by matching the ionisation rate for a thermal 
sample of He* to the results reported in Ref. [16] . 

The important dynamics of this system occur inside or 
near the condensate, and the evolution far below this re- 
gion is well described by free-fall. We therefore employ a 
two-step method where we first calculate the evolution of 
the three magnetic substates in a restricted region enclos- 
ing the condensate and then propagate the momentum 



flux density at the edge of the simulation region clas- 
sically to determine the profile on the detector located 
4 cm below the BEC. To take into account the estimated 
expansion due to mean-field repulsion in the beam itself 
we use an approximate method where we convolve the 
detector profiles with a Gaussian of width of 200 pm. 

The results of these simulations are shown in the lower 
row of Fig. [I] and agree well with the experimental re- 
sults in the middle row. As borne out by these simu- 
lations, the process forming the peaks is the scattering 
of stationary trapped and untrapped atoms, which pro- 
duce pairs of atoms with oppositely-directed momenta 
in the weak trapping dimension. This is different to the 
usual FWM process where two atoms with nonzero ki- 
netic energy collide, conserving kinetic energy to scatter 
into different momentum modes. 

To determine the source of the additional kinetic en- 
ergy, we consider that the scattered atoms will cause in- 
terference with any background atoms. This interference 
causes density gratings for both of the internal states. In 
the case of equal scattering lengths the interaction energy 
is unchanged, as the total density is unchanged. However, 
He* has a reduced scattering length for 0-0 collisions as 
compared to 0-1 and 1-1 collisions, so the interaction 
energy is reduced by the formation of the gratings. In 
the case of equal, uniform background densities in both 
states, the change in the interaction energy is 

A(^i n t> =U(1 -k) J (2\^\ 2 ecos(kx) - e 2 cos 2 (kx)^j dx, 

(3) 

where |\I/| 2 is the background density of both atomic 
fields, e is the amplitude of the density gratings in both 
states, U = Un = Uio an d k = «o,o/^i,i < 1 is the ratio 
of the scattering lengths such that kU = [Too- The first 
term in this expression averages to zero over one period 
of the density grating and the second term is always neg- 
ative, indicating that the interaction energy has reduced 
due to the scattering. Thus the increased kinetic energy 
after scattering is balanced by the reduced mean field en- 
ergy, allowing the scattering to take place, even though 
the atoms were initially at rest. A resonance between 
this reduced interaction energy and the increased kinetic 
energy of the scattered atoms causes well-defined beams 
to be produced by FWM. 

Generalising this argument to allow different density 
backgrounds for the m = 1 and m = states yields 
the wavenumber corresponding to the energy-momentum 
resonance of the scattered atoms as 

k=^2MU (2|* 1 ||* | " |*i| 2 -«|*o| 2 )- (4) 

Figure [3] shows the result of a simulation of the density 
of the BEC after a period of outcoupling showing the 
appearance of this density grating in the trapped BEC. 
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FIG. 3: (color online) A Gross-Pitaevskii simulation of the 
in-trap condensate density after 3.6 ms of out coupling, mag- 
nified to show the central region of the condensate, where the 
scattering-induced grating has high visibility. The weak and 
tight trap directions are x and y respectively. 



The FWM process discussed in this Letter can be ei- 
ther stimulated by the presence of an existing atomic 
population in the final momentum states, or seeded spon- 
taneously if those target states are initially empty. 

The stimulated process occurs if the outcoupling sur- 
face is chosen to be an ellipsoidal shell, for example by 
detuning the outcoupling from the center of the con- 
densate by 2 kHz. This causes a range of densities in 
the trapped and untrapped fields to exist over the sur- 
face. From Eq. Q we see that this corresponds to a 
range of momenta that are potentially resonant, from 
zero up to some maximum value. Atoms accelerated in 
the mean-field and trap potentials are amplified by this 
FWM scattering, allowing resonant stimulated scatter- 
ing into higher and higher final momentum modes. This 
sweeps out the halo shown in Fig. [I] and finally reaches 
the momentum corresponding to the peaks which then 
become heavily populated. 

Spontaneous processes cannot be described by a mean- 
field model, so we simulated this system using the 
stochastic truncated Wigner (TW) method [18] . As the 
detuning is reduced below 2 kHz, the outcoupling surface 
decreases in size and the minimum resonant momentum 
hk for the FWM process increases while the maximum 
resonant momentum remains approximately constant. In 
this situation the lower momentum modes near zero are 
no longer resonant, and the FWM process can only be 
started spontaneously. After the spontaneous scattering 
has taken place, these momentum modes then build up 
over time through stimulated scattering. At these lower 
detunings, in the TW simulations the extra peaks are 
better defined due to the narrower range of resonant mo- 
menta (though weaker and more difficult to see experi- 
mentally). Importantly, the peaks are non-existent us- 
ing semiclassical simulations, where spontaneous scatter- 
ing does not occur. Moreover, the resonant momentum 
for the gratings of 1.2 x 10 6 m _1 is in excellent agree- 



ment with the result given by Eq. (jl} of 1.3 x 10 m . 
This demonstrates that the momentum modes are gen- 
erated by a spontaneously-seeded four-wave mixing pro- 
cess, where two atoms in the condensate with near-zero 
momentum are scattered into modes with opposite mo- 
mentum along the weak trapping axis. This process will 
produce EPR entanglement between the trapped and un- 
trapped scattered atoms. This may lead to entanglement 
between atoms on opposing points of the cones in Fig. [2j 
but due to the output dynamics and the vertical integra- 
tion of the detection scheme, this would be difficult to 
measure. However, number difference squeezing between 
the opposing peaks should be far more robust. 

As the maximum wavenumber for these resonances in 
both the spontaneously and coherently seeded process 
corresponds to a wavelength larger than the Thomas- 
Fermi radius in the tight trapping dimension, the en- 
ergy and momentum resonances only co-exist in the weak 
trapping dimension, leading to the asymmetry in Fig. [2] 
If the density in the untrapped state is too low, no grating 
will form as the wavenumber for the energy-momentum 
resonance given by Eq. Q will be imaginary as is the 
case for low outcoupling (upper row Fig. [T]) . 

We have shown that appropriate outcoupling from a 
He* BEC can produce well-defined additional peaks in 
the output beam. Field theoretic and semiclassical mod- 
els show that these peaks are formed from scattering of 
pairs of atoms in BEC, and are therefore entangled upon 
formation. It remains to be seen whether useful entan- 
glement remains after the outcoupling process. These 
beams are produced directly via RF outcoupling in a 
magnetic trap. The potential advantages of these cor- 
related beams are that they are spatially well separated 
from the background of the atom laser and that the quasi- 
continuous dilute beam will likely remain coherent over 
larger timescales than trapped fields. 
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